Europhysics Letters 



PREPRINT 



Boundary effect on CDW: Friedel oscillations, STM image 

Balazs Dora 

Department of Physics, Budapest University of Technology and Economics, 
H-1521 Budapest, Hungary 



PACS. 71.45.Lr - Charge-density-wave systems. 
PACS. 68.37.Ef - Scanning tunneling microscopy. 
PACS. 73.40.Gk - Tunneling. 



Abstract. - We study the effect of open boundary condition on charge density waves (CDW). 
The electron density oscillates rapidly close to the boundary, and additional non-oscillating 
terms (~ hi(r)) appear. The Friedel oscillations survive beyond the CDW coherence length 
(vf/A), but their amplitude gets heavily suppressed. The scanning tunneling microscopy image 
(STM) of CDW shows clear features of the boundary. The local tunneling conductance becomes 
asymmetric with respect to the Fermi energy, and considerable amount of spectral weight is 
transferred to the lower gap edge. Also it exhibits additional zeros reflecting the influence of 
the boundary. 



Introduction. The behaviour of charge density waves (CDW) in bulk systems has 
received considerable attention in the last few decades [1]. But unlike mesoscopic supercon- 
ductors, only a few attempts have been made to explore the response of CDW in mesoscopic 
systems [2-5]. These mainly concern charge density wave junctions, and concentrate on tun- 
neling between CDW and other systems. In this letter, we are going to study the effect of an 
open end along the CDW chain on the electron density and scanning tunneling microscope 
(STM) image, as shown in fig. ^ 

As a model we consider a quasi-one dimensional charge density wave, treated in the mean- 
field approximation. Due to quasi-one dimensionality, the effective mass of the electrons in 
the interchain direction is much larger than along the chains, and it is reasonable to assume 
that the main affect of neighbouring chains is to suppress the thermal fluctuations of the order 
parameter [1,4]. Impurities in low dimensional systems tend to cut the sample. So instead of 
introducing an external potential representing the impurity, we shall mimic its effect by a semi 
infinite chain, namely by open boundary condition [6]. Hence the present model is thought 
to describe impurities in the strong pinning limit, or the behaviour of finite CDW with open 
boundaries. Similar models in one dimensional systems have extensively been studied over 
the years [7-9]. 

Experimentally, in-chain tunneling studies have been made on CDW systems [10, 11], and 
most of STM scans were performed on bulk CDW materials such as Fe doped NbSe3 [12] or 
NbSe2 [13], where boundary condition or single impurity results were not reported. CDW 
film growth and structuring have been started from oxide Rbo.3oMo03 (blue bronze) [11, 14] 
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and NbSe3, o-TaS3 [15], but no exhaustive studies of STM imaging were made so far close to 
interfaces or insulating barriers. Friedel oscillations in vanadium-doped blue bronze [16] have 
successfully been detected around V substituant. 

In a normal metal with open boundary (semi infinite chain) , the electron density behaves 

as 



where r is measured from the open end. It vanishes right at the boundary, and produces 
density oscillation with a periodicity of 7r/fcp. In the followings we are going to study what 
happens to a CDW close to an open boundary. 

Friedel oscillations. - For simplicity we consider a system of spinless fermions in CDW 
state. The inclusion of spin does not alter our results. The effect of open boundary can 
readily be incorporated into the theory by requiring, that the field operators should vanish at 
the origin: a r= o = 0- This can be fulfilled by connecting the states with different wavenumber 
but equal energy [7], as sketched in fig. [21 Mathematically, this means, that in the reduced 



Fig. 1 - Schematic view of the experimental setup. The STM tip moves along the CDW sample. 

Fig. 2 - The schematic energy spectrum of CDW is plotted. Due to the boundary condition, left and 
right movers are connected as shown by the red dashed line. 

Brillouin zone (0 < k < 2kp), the new quasiparticle operators, which diagonalize the CDW 
system [1,17], fulfill the relations 




(1) 




Fig. 1 
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(2) 
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and the original electron operators are expressed as 



a-k-Q 



e^iukd+.k + UfccL, fc) 
Vkd +i k - Ufcrf-,fe 



(4) 
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with 

u * \ l (l± m \ (5) 

Q — 2kp is the nesting vector, kF is the Fermi wavenumber, £(fc) = Vp(k — hp), vf is the 
Fermi velocity, A is the CDW order parameter. The phase of the CDW is locked to its 
optimum value determined by the proximity of the boundary, namely cj> — 0. This assures 
that the CDW contribution to the electron density at the boundary is minimal. The very 
same results can be read off from refs. [18,19], where the effect of point-like impurities were 
explored. Since we are not interested in phenomenon involving collective modes such as the 
sliding of CDW, only static (electron density) quantity and tunneling perpendicular to the 
chain will be considered, the above simplification is well motivated. Due to the proximity 
of boundary, A(x) is not constant close to the boundary, but decays smoothly over a finite 
distance from the boundary. However, the phenomena discussed here are not expected to be 
altered significantly by position dependent corrections to A [4,5]. 
The electron density is obtained as 



n(r)=2 £ W(Q&) 

0<fc<2fc- L \ f / 



1 1 cos 2 (k F r) + sm 2 (k F r) cos ' 



Vf 



_sw_ sin ( 2 fc F r) - A cos(2k F r) - cos sv ' 



— tanh £(k) sin — sin(2fc F r) — A cos(2/c^r) — cos 

2E 2 \ \ vf J V V v f 



(G) 



where r is measured from the open boundary, E — y/£,(k) 2 + A 2 . With the use of the gap 
equation [17], this can be further simplified. At T = 0, the total density is evaluated as 

k F 2A , , A 

n{r) = — cos(2fc F r) + F(2k F r), (7) 

7r g ttvf 

where 

F(x) = K (0\ + Ci(x) - sin^Kx (^j , (8) 

Ci(x) is the cosine integral, K n (x) is the nth Bessel function of the second kind, W = vpkp, 
g > is the effective electron-phonon coupling. The first term in eq. is the homogeneous 
electron density, the second one describes the spatially periodic charge density oscillations, 
while the F(x) function contains all the information concerning the effect of the boundary. 
It is shown in figs. |I]for A/W = 0.01 and 0.1 together with the oscillations in a normal 
metal. ni, oun dary / 'tiq = AF(2kpr)/W , uq = kp/ir is the density in a homogeneous system. 
The difference between a normal metal and CDW becomes more pronounced for larger values 
of A/W. 

Right at the boundary, n(r = 0) = as required. In the l/kp <C r <C vf/A region, the 
Friedel oscillations caused by the boundary decay as 

n A A v f \ sin(2fc F r) 

nboundary(X) = I hi — T 7 ) , (9) 

ttvf v rA ) 2nr 



where 7 = 0.57721 is the Euler's constant. The first term represents the CDW contribution, 
while the second one is identical to that in a normal metal [18] as shown in eq. ^ Spatially 
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Fig. 3 - The Friedel oscillations caused by the boundary are shown in a CDW for A/W = 0.01 (solid 
line) and in a normal metal (dashed-dotted line), no = fcj? /tt is the density in a homogeneous system. 
The red dashed line and green circles denote the asymptotic formulas for small and large r. 
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Fig. 4 - The Friedel oscillations caused by the boundary are shown in a CDW for A/W = 0.1 (solid 
line) and in a normal metal (dashed line). Due to the shorter coherence length, the oscillations die 
out at shorter distances as opposed to fig. 



oscillating terms arise also from CDW, but their amplitude are typically A/W times smaller 
than the second term in eq. 
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On the other hand, beyond the CDW coherence length, 

A sin(2fci?r) 

^boundary y ) — • l-^J 

TTV F 2fei?r 

This means, that the amplitude of the Friedel oscillations changes by a factor of A/W as one 
passes through the CDW coherence length from the boundary. This makes the detection of 
such fine structures superimposed on the usual charge density oscillations very difficult. On 
the other hand, the short distance behaviour could readily be checked experimentally. A more 
direct way to detect boundary effects is to look for probes with energy and spatial resolution 
such as the local tunneling conductance, as will be discussed in the followings. 




HLmmmm* 

— - 
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Fig. 5 - The local tunneling conductance is shown close to the boundary as a function of voltage (V) 
and distance from the boundary (r) for W/A = 20. The tunneling current along the horizontal black 
line is plotted in fig. |S| 

STM image, tunneling conductance. - The STM current I is directly related to how many 
electrons are locally available in the CDW. At a position r, the local tunneling conductance 
measures the local density of states, and is given as 



dl(r) 
dV 



N(oj,r) 



1 

2^ 



e luJt ({a r (t),a+(0)})dt. 



(11) 



Usually the tunneling matrix element between the sample and the tip depends on the wavevec- 
tors of the electron in the CDW and in the tip. However, the behaviour of CDW is mainly 
determined by electrons living in the Fermi surface, which is determined by the wavevector 
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component corresponding to the quasi-one dimensional direction. In this sense, this com- 
ponent carries all the informations about the condensate, and the inclusion of perpendicular 
components is not expected to alter the results. After straightforward calculation ea.llllvields 
to 



dl(r) 4|w| 



dV iry/uj 2 ~ A 



= ^sin 2 (\/aj 2 — A 2 r/vp) cos 2 (/cf r) + sin 2 (kpr) cos 2 (\/ lu 2 — A 2 r/vp)+ 



Vlu 2 - A 2 . 



H sin(2\/ lu 2 — A 2 r/vp) sm(2kpr) H (cos(2kpr) — cos(2\/u; 2 — A 2 r/vp)) 

2lo 2lo 



(12) 



with E = \J (J 1 — A 2 . In a normal metal this reduces to 



dl(r) 4 2 f tor 



— sin 



— + k F r). (13) 



dV 7T \Vp 

Hence the position of zeros in both cases is determined as r — Vprm/{u) + W), n a natural 

number. As u> increases, the pattern gets denser, as can be checked in fig. [5] Along the 

black line, a typical plot of the local density of states and the tunneling current is shown in 

fig. El The tunneling current remains unchanged with varying voltages at the zeros of the 

vpfiir 

local density of states, namely at ui = W . In a homogeneous CDW, we expect sharp 

r 

peaks at both A and —A. In the present case, however, due to the presence of bound states 
caused by the boundary, only one peak remains situated at —A, in accordance with ref. [19]. 
Here the effect of a single impurity was studied. After letting the strength of impurity to go 
to infinity, the bound state induced by the scatterer moves to —A, and the divergent peak at 
A disappears. 

Conclusion. We have investigated the effect of an open boundary on charge den- 
sity waves. It is believed that the effect of strong impurities or tunneling barriers can be 
approximated by an open end. The right and left moving particles mix up. The electron den- 
sity consists of the usual CDW contribution and terms caused by the proximity of boundary. 
Within the CDW coherence length (~ vpj A), these terms correspond to the response of a free 
electron gas and additional ln(r) corrections caused by the interaction between the boundary 
and the condensate. Beyond the CDW coherence length, only the small amplitude oscillations 
survive, but their detection seems to be a very difficult task to deal with. Quantities with 
both energy and spatial resolution might better help to clearly see boundary effects. Among 
them we have chosen to study the local tunneling conductance along the chain, measurable 
by scanning tunneling microscope. The differential conductance is zero for voltages smaller 
than the gap maximum. For negative V, dI(r)/dV exhibits a sharp peak ~ Aj \JV 2 — A 2 as 
V approaches —A. For positive voltages, the differential conductance increases smoothly with 
V, and no divergences are found. This is in accordance with ref. [19]. The tunneling current 
changes inflection with increasing tip voltage due to the presence of zeros in the local density 

of states of CDW. Also it remains unchanged not only around the Fermi energy, but also at 

Vpn-K 

the additional zeros of the local density of states given by u> = W. 

r 
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Fig. 6 - The STM tunneling current (I(r)) is shown at kpr = 8 as a function of the applied voltage 
along the black line in fig- EI The inset shows the tunneling conductance (dI(r)/dV) along the same 
line. 
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